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Abstract. 

We derive exact relations for SU(2) lattice gauge theory in 3+1 dimen- 
sions. In terms of Abelian projection, these are the expectation values 
of Maxwell equations that define a new field strength operator and 
conserved, dynamic electric currents formed from the charged matter 
and ghost fields. The effect of gauge fixing is calculated, and in the 
maximally Abelian gauge we find antiscreening of U(l) Wilson loop 
source charges. We discuss the importance of these quantities in the 
dual superconducting vacuum mechanism of confinement. 
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1 Introduction 



Lattice studies based on Abelian projection have had considerable success identifying the 
dynamical variables relevant to the physics of quark confinement. There is no definitive 
way as yet of choosing the optimum variables, but in the maximally Abelian gauge [[[], 3 
the U(l) fields remaining after Abelian projection produce a heavy quark potential that 
continues to rise linearly [Q]. Further the string tension is almost, but not exactly, equal to 
the full SU(2) quantity; 92% in a recent study at (3 = 2.5115 g. 

This suggests that we may be close to identifying an underlying principle governing 
confinement. All elements of a dual superconducting vacuum appear to be present |5|, 
in the maximally Abelian gauge magnetic monopoles reproduce nearly all of the U(l) 
string tension |], ffl. The spontaneous breaking to the U(l) gauge symmetry is signalled 
by the non-zero vacuum expectation value of monopole operator [[?], [|. The profile of the 
electric field and the persistent magnetic monopole currents in the vortex between quark and 
antiquark are well described by an effective theory, the Ginzburg-Landau, or equivalently 
a Higgs theory giving a London penetration depth and Ginzburg-Landau coherence length 



1C 



Central to finding the effective theory is the definition of the field strength operator in 
the Abelian projected theory, entering not only in the vortex profiles but also in the formula 
for the monopole operator. All definitions should be equivalent in the continuum limit, but 
use of the appropriate lattice expression should lead to a minimisation of discretisation 
errors. 

In this paper we exploit lattice symmetries to derive such an operator that satisfies 
Ehrenfest relations; Maxwell's equations for ensemble averages irrespective of lattice arte- 
facts. In section ^| we introduce and review this method in pure U(l) theories [flUI , dis- 
cussing the Abelian projected SU(2) theory, with and without gauge fixing to the maximally 
Abelian gauge, in section ||. 

The charged coset fields are normally discarded in Abelian projection, as are the ghost 
fields arising from the gauge fixing procedure. Since the remainder of the SU(2) infrared 
physics must arise from these, an understanding of their role is central to completing the 
picture of full SU(2) confinement. In section ^| we begin to address this issue, showing 
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that these fields form a charged U(l) current, and in section f| demonstrate that in the 
maximally Abelian gauge the supposedly unit charged Abelian Wilson loop has an upward 
renormalisation of charge due to this current of 0(15%). A localised cloud of like polarity 
charge is induced in the vacuum in the vicinity of a source, producing an effect reminiscent 
of the antiscreening of charge in QCD. In other gauges studied, the analogous current is 
weaker, and acts to screen the source. 

We show that this current can be quantitatively written as a sum of terms from the coset 
and ghost fields. The contribution of the ghost fields in the maximally Abelian gauge in 
this context is found to be small. The effect of the the Gribov ambiguity on these currents 
is argued to be slight. 

Finally in section |5| we discuss these results. Some preliminary results have already 
appeared |12|. The reader's attention is drawn to related work in this subject |13], [14|. 



2 Abelian theories 



The Wilson action 



Sw = ( 1 ~~ cos 6^(71)) 



(1) 



comprises link angles {0 M (n) G [— tt, 7t)} summed to form plaquette angles 9 i 



ft]' 



n 



0» + 

9 v (n+ji) — M (n+z>) — 9 v (n). External electric sources may be represented by a Wilson loop, 
and we consider here for simplicity a specific plaquette on the lattice P K \(n') = expi9 K \(n') 
with real and imaginary parts R K \, I K \ respectively. The partition function 



J SKC 



[de„]P KX (n')e-^. 



(2) 



is invariant under the introduction of an arbitrary constant into any link angle. We call 
this a shift invariance and we focus on the consequence of this on a particular link: M (n) — > 
@ii{ n ) + e - The shift corresponds to either a left or right multiplication of that link by a 
group element: e^^™* 1 — > e^^e* 6 = e t£ e t9 ^( n \ The Haar measure is invariant under this 
transformation: d(8^(n) + e) = d6^(n). 

The first order shift in the partition function gives an identity: 



1 dZ. 



snc 



J S11C 



de 



e=0 



■'STIC 



dP K \(n') 
de 

2 



dS 



£=0 



de 



e=0. 
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f[d9r]R KX .e-f )S * 

if3 



iSw — 



(Rkx) 



dS w 



iR K \.5 w - H K \.(3 



de 



e=0. 



-PS W 



de 



(3) 



£ = 0/ 



where 5w 7^ only when the shifted link coincides with a link in the source, 



(4) 



The real or imaginary part of P K \ is dropped if it contributes a term odd in the link angle 
and hence has a zero expectation value. 

Multiplying eqn. (|J) by the electric charge e, where (3 = 1/e 2 , we use the backwards 
lattice difference operator to define the lattice field strength tensor, f^ v 



i as, 



ir 



e de 



e=0 



The identity becomes: 



(I KX .A-f^(n)) 



(Rkx) 



where the static current density is 



j^(n) = e5 w 



(5) 



(6) 



(7) 



We have arrived at what appear to be a discretised version of the continuum Maxwell's 
equations for the U(l) fields, but satisfied by the expectation values rather than merely in 
the classical limit of extremising the action. In quantum mechanics, Ehrenfest's theorem 
relates the time derivative of the position operator to the potential in a way reminiscent of 
Newton's classical equations of motion, (x) = — (VV(x)) . By analogy, we label the 
lattice expressions 'Ehrenfest identities'^ 

The theorem defines the lattice field strength operator whose form is dictated by the 
derivative of the U(l) lattice action. As all actions (in the same universality class) are 
equivalent in the weak coupling limit, so too are the corresponding definitions of the field 
strength operator. Much work, however, is performed at finite lattice spacing, and it is 

1 A term also used in Zach et. al. 
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advantageous to avoid extraneous 0(a 2 ) effects by using the 'correct' operator. Using this 
operator we may measure the charge density, which the theorem shows is exactly that of 
the introduced source. No further charge is induced. 

3 Non— Abelian gauge theories 

The correct Abelian field strength operator could be calculated as above if the effective 
Abelian action for the U(l) fields after Abelian projection were known. It is not, and 
we approach the problem from the full SU(2) action. The SU(2) theory has symmetries 
analogous to those of the pure Abelian theory, each of which gives rise to Ehrenfest iden- 
tities. Since these, in the continuum limit, resemble the Euler-Lagrange equations for the 
corresponding continuum action, we begin by briefly considering these. 

The continuum Lagrangian is C = \G^ lu (x)G" xl/ (x) , the isospin index a e {1,2,3}. 
Under Abelian projection the third component of the gauge field becomes the Abelian 
gauge potential. We can rewrite £ to emphasis this: 

c = \ [f, u f, v + w; v w, v + 1 -f, u {w,w; - w;w u ) - 1 (w,w; - w;w u ) 2 ^j (8) 

where F lJiV {x) = d^A z u {x) —d v A? {x) is the Abelian field strength, the remaining components 
forming a complex matter field W^{x) = A^x) + iA 2 {pc). This field is electrically charged 
with respect to the photon; D fM W u (x) = (<9 M — iA* (x))W v (x), giving W fll/ (x) = D IM W L ,(x) — 
D u W,(x). 

Consider the extremisation with respect to A*{x). This gives what appear to be 
Maxwell's equations with a dynamical, real-valued, conserved electric current formed from 
the coset (matter) fields and their coupling to the photon. 

d„F, v {x) = j^(x) = ^ (2 (w, K w: - w; K w R ) - d K (w K w; - W* K W,)) . (9) 

The first two terms in the current are precisely what would be expected for a charged 
vector field, i.e. where D^W^x) = 0. This is not in general true. It is interesting to note, 
however, that the imposition of this constraint amounts precisely to fixing the theory to 
the maximally Abelian gauge. 
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The derivation here relied upon the assumption that we might vary the photon field 
independently of the charged coset fields. While this is true in the full theory, a gauge fixing 
constraint couples the variations in the fields. Such an extremisation problem is usually 
tackled using Lagrange multipliers. The correct lattice operators cannot be predicted by 
naive discretisation of continuum results, so we do not pursue this approach here but move 
on to the lattice Ehrenfest identities. 



3.1 Abelian projection on the lattice 

After gauge fixing, the SU(2) link matrices may be decomposed in a 'left coset' form: 



cos(0 M (n)) sin(0 M (n))e i > (n) \ 



V 



sin 



/> M (n))e- J >(«) cos(0») j 
Under a U(l) gauge transformation, {g(n) = exp [ia^n)^]}, 



e -™„ 



(n) 



(10) 



9/x(n) -> 9^{n) + a(n) - a(n + fi) ln{n) -> 7^( n ) + 2a(n) 



(11) 



In other words, the left coset field derived from the link U^n) is a doubly charged matter 
field living on the site n and is invariant under U(l) gauge transformations at neighbouring 
sites. 

The c M = cos(0 M ) are real-valued fields which near the continuum ~ 1 + 0(a 2 ) where 
a is the lattice spacing. The off-diagonal = sin(0 M )e nfl become the charged coset fields 
gaW^(x), and 6> M the photon field gaA^(x). [The SU(2) coupling f3 — in 3+1 dimensions.] 

The SU(2) shift symmetries are the left and right multiplications of a link by an arbitrary 
constant SU(2) matrix, under which the Haar measure is invariant. Since it is the a = 3 
component that becomes the photon, we consider here only shift matrices of the form 



right : U„(n) -> U„(n)U = {w^n) -> w^n) , 0„(n) -> 0„(n) + e} 
left : U^n) ^ UU^n) = {w^n) -> w^n)^ , 0» -> 0» + e} 



(12) 
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3.2 The identities — no gauge fixing 



We first derive the Ehrenfest relations in the simpler, but artificial, context of Abelian 
projection without gauge fixing, 'no gauge.' The SU(2) link matrices combine to form 
plaquettes U^in) = U^{n)U u {n + ji)UUn + u)Ul(n) and the Wilson action 

Sw= £ fi-^M"))' ( 13 ) 

Writing each link as the sum of a diagonal and an off-diagonal matrix, U = D + O, we 
define the product of the diagonal terms around a plaquette 

Qnx{n') = D K (n')D x (n' + k)D{(n' + X)D[(n'). (14) 

and choose as a source term a specific plaquette P K \(n') with real and imaginary parts: 

Rk\ = r Tr{Q KA (n')} = (c K (n').c A (n' + k).c K (n' + A).c A (n')) cos6» KA (n'), 
4a = 2 I^O^QreACn')} = (c K (ra').c A (ra' + k).c K (n' + A).c A (n')) sin6 l KA (n / ). (15) 

It is for consistency that we include the coset fields. The fields c K (n') — > 1 in the continuum 
limit and in the maximally Abelian gauge can be considered to be almost constant [ITR |To|, 



since the fluctuations are very small. For this reason, we anticipate only minor changes in 
the measured averages if the fields c M (n) were excluded from the source as is the case in 
the traditional Abelian source loop. 

The partition function is 

Zsnc = J [dU„]P KX {n')e-^ (16) 

If all fields bar the 6^ are discarded, the SU(2) action reduces to the U(l) Wilson action, so 
we define an effective electric charge, e, by (3 = \. The shift invariance gives an analogous 
result to the U(l) theory 



= (17) 

e=0/ 



(R.x) \ M de 

and the static current density is also based on the localised form of eqn. (^): j^ static ) (n) = 
e5w Although we have introduced an Abelian projected source, the links, U^, in the 
derivative of the action are SU(2) matrices. Were we to replace them also with their 



diagonal components only, = diag (c^e* 6 ^ , c M e ldfJ -), we should get a result analogous to 
the U(l) case: 



IDS 



e de 



-A;V(n) = A;/^(n) (1* 



e=0;t/^D 

The remaining terms involving both D and O constitute a U(l) gauge invariant dynamical 
current. 

= Xf„ u (n)-£*»\n) (19) 



e de 



£=0 



The terms in jj 1 yn) that survive into the continuum limit yield the continuum Euler- 
Lagrange result (eqn. @). The Ehrenfest identity is thus: 



(I KX .A-f, u (n)) 



(RkX/ 



(20) 



3.3 Gauge fixing 

We restrict the discussion here to the maximally Abelian gauge, defined as the local max- 
imisation by gauge transformations {gin) = exp[ia a (n)cr a ]} of the Morse functional 



R{U] = -Y,n{ia 3 .Up(n).i<T 3 .Ul(n)}. 



(21) 



n,fi 



Denoting the gauge transformed link variable as U^{n) = (yf(n)[/ At (n)5'^(n + /t), at a local 
maximum, assumed to occur at g — 1, the first derivative with respect to gauge transfor- 
mations of R[U 9 } will be zero: 

da r(P) a=0 

where X(p) = E At >o ( u »(p)- i(T 3- u t(p) + #J(P - A)-^3-^(p - A)) • The residual U(l) sym- 
metry is seen in F 3 being trivially zero. The gauge fixed partition function is thus 



2e 3rc Tr {ia c .X(n)} = 



(22) 



Z 9 Jnc = l[dU,}P K ,[U-n'].e-^ u \/\ FP [U] (l[HF r [U;p]) 

\p,r 



(23) 



and App is the Faddeev-Popov gauge fixing operator, a Jacobian to reweight the integration 
measure after the introduction of the constraint into the partition function: 



A FP [U]= ( [ldg]H5(F r iU°;p]) 







f- 


det M; s q 



. njrs dF r [U°;p} 

and pq s ~^wr 



(24) 



a=0 



is the second term in the Taylor expansion of R |[7|. The integral expression is gauge 



invariant by the invariance of the Haar measure. The determinant is gauge invariant by 
fiat; we must evaluate it 'on the constraint' (i.e. at F = 0). If F ^ we must first move 
the configuration along the gauge orbit by gauge transformations until it does. 

Applying a shift to a single link and differentiating with respect to the shift parameter, 
e, now raises problems, as the shift drives the previously gauge fixed configuration off the 
constraint. This manifests as derivatives of the constraint 5-functions, which we must 
avoid. One approach would be to shift not one link, but each of a Polyakov line of links 
extending around the lattice. By careful choice of the relative sizes of these shifts we may 
conspire to remove the offending 5-function terms using the constraint equations. This 
is at the expense of the shift as a local probe of charge density. Instead we accompany 
the shift of one specific link by an SU(2) gauge transformation over the whole lattice that 
is also linear in the shift parameter. The combined effect of the shift and 'corrective' 
gauge transformation is to move on a path in configuration space parameterised by e that 
remains on the trajectory satisfying the maximally Abelian gauge. There is potentially 
also a non-locality here, but we shall see that it is very limited in its extent. 

The corrective gauge transformation {g(q) = 1 + ier] s (q)a s } is calculated using the in- 
verse of the Faddeev-Popov matrix, and so 7/3 = 0: 



Vs(q) = -1^(M ) qp — 

p,r ut 



(25) 

e=0 



We now apply both the shift and corrective gauge transformation to the partition function 
in eqn. (p3|), and differentiate with respect to e. The action is gauge invariant, and thus 
only the shift has an effect and the field strength operator and current jj? yn) are as before 
(eqn. ©). 

The derivative of the source plaquette now makes two contributions; the shift gives the 
same static current as before. The second comes from the corrective gauge transformation, 
and forms a conserved current 7'ff auge ). 

Finally, there are contributions that arise from perturbing the links making up the 
Faddeev-Popov operator. Under this, the matrix M — > M + eN. As we have remarked, the 
Faddeev- Popov operator is specifically the determinant of M evaluated on the constraint. 
By introducing rj we have stayed on the constraint as we shifted the link, and N = A + B 
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has contributions from the shift, A, and from the corrective gauge transformation, B. The 
derivative forms a further conserved current: 



,-(FP) _ 



e dA FP [U £ 



A FP [U] de 



= e 



ddet(l + eM- 1 N) 



6 = 



de 



e=0 



= eTr {M^N} 



(26) 



Defining j£ otal) = j^ dyn) + j<f augc) + jf p) , we have as the final result a set of identities 



- "gf_ _ j< static )( n ) _|_ x 



gf 



{Rn\) gf 

where the gauge fixed expectation value is defined as 



(27) 



gf 



(0) gf = j[du,}.o.e-^.A FP [u] {uH Fr [U;p})) ■ 



(28) 



We conclude this section with two remarks. In simulation, a lattice average using SU(2) 
configurations calculated using standard Monte Carlo methods and then gauge fixed will 
include the Faddeev-Popov operator in the importance sampling of the measure and it 
need not be explicitly calculated. 

Secondly, when deriving the Faddeev-Popov matrix, there is some ambiguity as to 
whether one should start from the constraint 



lvl pq 



dF r [U g ;p] 



da s (q) 

or from the original functional 



d f dR[U 9 ] 
a=0 da s (q) \ da r (p) 



a=0; 



(29) 



Q = 



M rs — d2R[U9] 
pq da s (q)da r (p) 



(30) 



a=0 



Since we believe that M should be a symmetric matrix in general, the latter seems the 
most natural. In practise, we find that the differences between the two approaches are 
all multiples of the constraint F and hence when the gauge condition is satisfied the two 
definitions coincide. Similar ambiguities arise when the response of the Faddeev-Popov 
matrix to the corrective gauge transformation is considered, since for this we require the 
third term, L, in the Taylor expansion of R 



u,t 



(31) 
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and it is unclear whether we should begin from R, F or M in its derivation. Again we find 
that the most symmetric tensor is derived from R, but that all expressions are the same 
on the constraint. This we believe to be a general property. 

The ambiguities are more serious in the case of the contribution of the shift to the 
Faddeev-Popov matrix. The matrix A is derived from R by two differentiations with 
respect to gauge transformations, and one with respect to the shift. The order of these, 
and at what point variables are set to zero does appear to matter in this case, even when 
F = 0. The correct order of derivatives is as follows: 



d ( d ( dR[U°] 
pq da s (q) \de\ da r (p) 



a=0/ 



e=0/ 



(32) 

Q=0 



The inner two nested derivatives give the first order shift correction to the constraint, i.e., 
the second term in F + eF t = 0. The Faddeev-Popov operator is the matrix of derivatives of 
this about the constraint and is formed by the outermost nested derivative. The corrected 
Faddeev-Popov operator must be calculated about the shifted constraint. If, for example, 
we reverse the order of the outer two nested derivatives we would be calculating the shift 
of the lowest order Faddeev-Popov operator. This would give an incorrect contribution to 
the current. Numerical studies confirms these conclusions. 



4 Numerical investigation 

The terms in the Ehrenfest identities may be measured using Monte Carlo simulation. We 
begin with a careful test of the expressions, and for this we rewrite eqns. ( |2"0"| , |2"TD to minimise 
the combination of statistical errors: 



dSw 
de 



£ J-k^(n).(R KX ) = 

where the currents k differ from j by a factor of the charge. We further split fc|f p) into the 
(normally summed) contributions from the shift and corrective gauge transformation. To 
measure the gauge fixed currents separately requires the inversion of the large Faddeev- 
Popov matrix. This limits us to the unphysical lattice size of 4 4 for this test, where although 
the identities must still hold exactly, there may be significant finite volume effects on the 
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individual currents. For comparison we also test the 'no gauge' expression on the same 
lattice. 

In Table [l] we show results where the (left) shifted link is part of the source plaquette 
(i.e. | ^(static) | _ -q rj^g identities are verified within very small statistical errors. The 
numbers are for right shifts are identical within these errors. Off the source, the identities 
are equally well satisfied. 

Table |2| shows the normalised currents arising from shifting a timelike link included 
in the source, and then one spacelike lattice spacing away from the source in the same 
plane. In the maximally Abelian gauge we see that on the source the charge of the Abelian 
projected Wilson loop has been renormalised from unity to a value 0(15%) higher. Further, 
in the vicinity of the source we find a cloud of like charges has been induced in the vacuum. 
This is reminiscent of the charge antiscreening (or asymptotic freedom) of the full gauge 
theory. The charge cloud is localised, falling to near zero by two lattice spacings from the 
source. We note that in the case of 'no gauge' (and also in the Polyakov and diagonal 
plaquette gauges e.g. 'F^) there is also a renormalisation, but that it reduces the charge 
of the Abelian Wilson loop. The induced charge cloud is weaker and acts to screen the 
source, at odds with the behaviour of the full non- Abelian theory. In U(l) we might expect 
renormalisation and screening, but not in the pure gauge theory. 

Assuming the Ehrenfest identities to be true, we may infer k^° tal) from a measurement 
of the derived field strength operator, and break it down into some of its components. This 
allows the consideration of lattices of a more physically interesting size. In Table ^| we 
see that the currents are remarkably stable and free of finite volume effects even moving 
to lattices large enough to support infrared physics. Finite volume effects on the currents 
induced by a plaquette-sized source would thus appear to be slight. 

The Gribov ambiguity was neglected throughout. The Morse functional R[U] typically 
has a number of local maxima, each giving a different set of maximally Abelian gauge 
fields corresponding to the Gribov copies. There is some ambiguity regarding which of 
these to use and one is usually selected randomly during gauge fixing. The variation of 
gauge variant observables, such as the string tension after Abelian projection, between these 
copies, although not zero, is small enough that we may continue to neglect the ambiguity in 



studying the properties of this gauge [18|. It is possible, however, that observables derived 
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directly from the Faddeev-Popov operator may be unduly sensitive to this ambiguity. 
Only on lattices comparable in size to the confining length scale (~ 1 fm) do we see 
nonperturbative effects, and it is on these lattices that R[U] begins to have multiple maxima 



and Gribov copies appear [ 18[] . The Ehrenfest identities were derived by infinitessimally (in 
principle) shifting a configuration away from a maximum of R[U], and applying a corrective 
gauge transformation to return to the maximally Abelian gauge. If the Gribov copies are 
not too numerous, it seems likely that the configuration does not forsake the attractive 
influence of the original maximum, and thus returns to the 'same' Gribov copy. When the 
copies become very numerous, however, it may be that two maxima of R are sufficiently 
proximate that there is a 'flat direction' between them, characterised by an extra (near) 
zero eigenvalue in the Faddeev-Popov matrix. 

A 4 lattice at (3 = 2.3 is far too small to support nonperturbative physics, and indeed 
the configurations we studied there exhibited only one maximum of R in 1000 gauge fixings 
(we differentiated the maxima using the value of R and the U(l) plaquette action after 
Abelian projection, as in |18[). In none of our simulations were there any problems inverting 
the Faddeev-Popov matrix, and we believe the effects of the Gribov ambiguity to be very 
small. Only below (3 ~ 2.15 do Gribov copies appear, and although we again had no 
apparent problems inverting, the statistical noise prevents us from making any statements 
about the currents. On larger lattices, we have shown in Table |3| that the finite volume 
effects are slight, even on moving to physically large lattices. This does not prove that the 
Gribov copies have no effect (since we only undertook one gauge fixing per configuration), 
but the continuing smallness of the statistical errors is perhaps indicative that that. 

In a perturbative treatment of the gauge fixed action, the non-local Faddeev-Popov 
operator would be replaced by a Gaussian integral over propagating ghost fields. In this 
sense, then, the Faddeev-Popov current that we have isolated may be regarded as the 
contribution of the ghost fields, hitherto little studied in the maximally Abelian gauge. We 
conclude that in this context at least these fields play little role since the Faddeev-Popov 
current is only a fraction of the other currents. 
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5 Summary 



In this paper we have exploited symmetries of the lattice partition function to derive a 
set of exact, non-Abelian identities which define the Abelian field strength operator and 
a conserved electric current arising from the coset fields traditionally discarded in Abelian 
projection. The current has contributions from the action, the gauge fixing condition and 
the Faddeev-Popov operator. Numerical studies on small lattices verified the identity to 
within errors of a few per cent. We have found the Faddeev-Popov current in particular 
to be unusually sensitive to systematic effects such as low numerical precision and poor 
random number generators, but the origin of any remaining, subtle biases, if they exist, is 
not clear; we have already considered all terms in the partition function. 

In a pure U(l) theory the static quark potential may be measured using Wilson loops 
that correspond to unit charges moving in closed loops, as demonstrated by | {A~ f Ufl ) \ = 8 W . 
In Abelian projected SU(2) the same measurements in the maximally Abelian gauge yield 
an asymptotic area law decay and a string tension that is only slightly less than the full 
non-Abelian value. In other gauges it is not clear that an area law exists — certainly it is 
more troublesome to identify. 

We have seen that in the context of the full theory the Abelian Wilson loop must be rein- 
terpreted. The coset fields renormalise the charge of the loop as measured by \(A~ f u ^)\ and 
charge is also induced in the surrounding vacuum. Full SU(2) has antiscreening/asymptotic 
freedom of colour charge, and in the maximally Abelian gauge alone have we seen analo- 
gous behaviour, in that the source charge is increased and induces charge of like polarity 
in the neighbouring vacuum. Whether this renormalisation of charge can account for the 
reduction of the string tension upon Abelian projection in this gauge is not clear. In other 
gauges, where Abelian dominance of the string tension is not seen, the coset fields appear 
to have a qualitatively different behaviour, acting to suppress and screen the source charge. 

In conclusion, the improved field strength expression defined by the Ehrenfest iden- 



tity does not coincide with the lattice version of |L9| of 't Hooft's proposed field strength 
operator [2t|. The Abelian and monopole dominance of the string tension invites a dual 
superconductor hypothesis for confinement. If this is to be demonstrated quantitatively 
such as by verification of a (dual) London equation then a a careful understanding of the 
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field strength operator is required. The Ehrenfest identities may provide this 



in 
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ON THE SOURCE 


P 


= 2.3 


P 


= 2.5 






'no gauge' 


MA gauge 


'no gauge' 


MA gauge 


# meas. 




10000 


10000 


10000 


20000 






0.0765 (3) 


0.6162 (8) 


0.0815 (2) 


0.68275 (58) 


fc (statlc)_ 


) 


0.0762 (1) 


0.5607 (5) 


0.0818 (1) 


0.63069 (21) 


(I KX .k^ C >) 






0.0469 (1) 




0.04463 (5) 




- shift 




0.0063 (1) 




0.00565 (3) 


(Lx-k^) - 


- corr. 


g.t. 


0.0034 (11) 




0.00133 (51) 


LHS eqn. ([[ 


3 


0.0003 (4) 


-0.0011 (17) 


-0.0003 (3) 


0.00045 (71) 



Table 1: Numerical tests of the Ehrenfest identities on 4 4 lattices, with and without gauge 
fixing. 



ON THE SOURCE 


P = 


2.3 


P = 


2.5 




'no gauge' 


MA gauge 


'no gauge' 


MA gauge 


^(dyn) 


-0.1266 (40) 


0.0800 (5) 


-0.1293 (25) 


0.0781 (7) 


h (gauge) 




0.0836 (2) 




0.0711 (1) 


fcCFP) 




0.0173 (20) 




0.0112 (3) 


(static) _|_ 1^ (total) 


0.8734 (40) 


1.1809 (20) 


0.8707 (25) 


1.1604 (10) 



OFF THE SOURCE (3 = 2.3 /3 = 2.5 

'no gauge' MA gauge 'no gauge' MA gauge 

fc (dyn) -0.0074 (40) 0.0139 (5) -0.0126 (40) 0.0145 (5) 

A; (gaugo) 0.0042 (1) 0.0037 (1) 

k (FP) 0.0001 (25) 0.0012 (3) 

^(static) + fc( totai) -0.0074 (40) 0.0182 (25) -0.0126 (40) 0.0194 (10) 

Table 2: Ehrenfest currents {I K \.k) / (R K \) on 4 4 lattices, with and without gauge fixing 
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ON 


THE SOURCE 


J (total) 


j'(dyn) 


j (gauge) _|_ j'( FP ) 


P = 


2.3, L = 


4 


0.1809 (20) 


0.0800 (5) 


0.1009 (25) 


P = 


2.3, L = 


6 


0.1776 (10) 


0.0674 (5) 


0.1102 (5) 


(3 = 


2.3, L = 


8 


0.1735 (39) 


0.0642 (25) 


0.1093 (41) 


P = 


2.3, L = 


10 


0.1840 (83) 


0.0654 (30) 


0.1186 (94) 


P = 


2.5, L = 


4 


0.1596 (7) 


0.0781 (7) 


0.0824 (13) 


P = 


2.5, L = 


6 


0.1610 (15) 


0.0757 (17) 


0.0853 (17) 


P = 


2.5, L = 


10 


0.1658 (40) 


0.0753 (16) 


0.0905 (35) 



Table 3: Finite volume effects on the Ehrenfest currents, {I K \.k) g ^ / (R K \) g j in the maxi- 
mally Abelian gauge. 
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